Decoherence and gate performance of coupled solid state qubits 
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Solid state quantum bits are promising candidates for the realization of a scalable quantum com- 
puter. However, they are usually strongly limited by decoherence due to the many extra degrees of 
freedom of a solid state system. We investigate a system of two solid state qubits that are coupled 



via cri*' (8> cri-" type of coupling. This kind of setup is typical for pseudospin solid-state quantum 
bits such as charge or flux systems. We evaluate decoherence properties and gate quality factors in 
the presence of a common and two uncorrelated baths coupling to a^, respectively. We show that 
at low temperatures, uncorrelated baths do degrade the gate quality more severely. In particular, 
we show that in the case of a common bath, optimum gate performance of a CPHASE gate can 
be reached at very low temperatures, because our type of coupling commutes with the coupling to 
the decoherence, which makes this type of coupling attractive as compared to previously studied 



proposals with ay 
CNOT gate. 



Jj) 



-coupling. Although less pronounced, this advantage also applies to the 
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I. INTRODUCTION 



Quantum computation has been shown to perform cer- 
tain tasks much faster than classical computers 0, D . 
Presently, very mature physical realizations of this idea 
originate in atomic physics, optics, and nuclear magnetic 
resonance. These systems are phase coherent in abun- 
dance, however, scaling up the existing few-qubits sys- 
tems is not straightforward. Solid state quantum com- 
puters have the potential advantage of being arbitrar- 
ily scalable to large systems of many qubits [|[ |[ D- 
Their most important drawback is the coupling to the 
many degrees of freedom of a solid state system. Even 
though recently there has been fast progress in improv- 
ing the decoherence properties of experimentally realized 
solid state quantum bits ^, ^ |ll), this remains a 
formidable task. 

Quite a lot is known about decoherence properties of 
single solid state qubits, see e.g. [Ol 13, 111], but much 



less is known about systems of two or more coupled 
qubits ||l^, |l^, |l^. However, only for systems of at least 
two qubits, the central issue of entanglement can be stud- 
ied. The physically available types of qubit coupling can 
be classified as Heisenberg-type exchange which is typ- 
ical for real spin-l/2-systems, and Ising-type coupling, 
which is characteristic for pseudospin setups, where the 
computational degrees of freedom are not real spins. In 
the latter, the different spin components typically cor- 
respond to distinct variables, such as charge and flux 
| pO| , 18 whose couplings can and have to be engineered 
on completely different footing. Previous work jl^, |l^ 

presented the properties of a system of two coupled solid 

(?) 

ay type cou- 



state qubits which are coupled via cry 
pling as proposed in H as the current-current coupling 



of superconducting charge quantum bits. 

On the other hand, many systems such as inductively 
coupled flux qubits |^ , capacitively coupled charge qubits 
0, 0] and other pseudo-spin systems ||l^ are described by 

a oi^ (g) ai''^ Ising-type coupling. This indicates that the 
computational basis states are coupled, which i.e. in the 
case of flux qubits are magnetic fluxes, whereas (J^/y are 
electric charges. The (T^-observable is a natural way of 
coupling, because it is typically easy to couple to. We 
will study a two qubit-system coupled this way that is 
exposed to Gaussian noise coupling to az, the "natural" 
observable. This e.g. accounts for the crucial effect of 
electromagnetic noise in superconducting qubits. We will 
compare both the case of noise that affects both qubits in 
a correlated way and the case of uncorrelated single-qubit 
errors. We determine the decoherence properties of the 
system by application of the well known Bloch-Redfield 
formalism and determine quality factors of a CNOT gate 
for both types of errors and feasible parameters of the 
system. 



II. MODEL HAMILTONIAN 

We model the Hamiltonian of a system of two qubits, 
coupled via Ising-type coupling. Each of the two qubits 
is a two-state system that is described in pseudo-spin 
notation by the single-qubit Hamiltonian pSl] 



1 
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where e is the energy bias and A the tunnel matrix ele- 
ment. The coupling between the qubits is determined by 



(1) 



.(2) 



*Electronic address: 3torcz@theorie.physik.um-muenchen.de 



an extra term in the Hamiltonian Hqq — 
that represents e.g. inductive interaction (directly or via 
flux transformer) in the case of flux qubits |, |0). Thus, 
the complete two-qubit Hamiltonian in the absence of a 
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dissipative environment reads 
1 w.-^ 1 



H 



2qb 



E 



(2) 



The dissipative (bosonic) environment is conveniently 
modeled as either a common bath or two distinct baths 
of harmonic oscillators, coupling to the cr^ components of 
the two qubits. This approach universally models baths 
which produce Gaussian fluctuations, such as the noise 
from linear electrical circuits. An example for a situ- 
ation described by a common bath is long correlation 
length electromagnetic noise from the experimental en- 
vironment or noise generated or picked up by coupling 
elements such as flux transformers Q. Short correlation 
length radiation or local readout and control electronics 
coupling to individual qubits p3| might be described as 
coupling to two uncorrelated baths of harmonic oscilla- 
tors. 

One should note that if the number of qubits is in- 
creased to more than two, there might also occur dissi- 
pative effects which neither affect all qubits nor only a 
single qubit but rather a cluster of qubits thus enhancing 
the complexity of our considerations 

In the case of two uncorrelated baths, the full Hamil- 
tonian reads 
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where each qubit couples to its own, distinct, harmonic 
oscillator bath H^. , i = 1,2, via the coupling term 
^Wj^(j)^ i — 1,2, that bilinearly couples a qubit to the 
collective bath coordinate X^*) = ^i^^i^- We again 
sum over the two qubits. In the case of two qubits cou- 
pling to one common bath, we model our two qubit sys- 
tem with the Hamiltonian 
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where H b denotes one common bath of harmonic oscil- 
lators. 

The appropriate starting point for our further anal- 
ysis is the singlet/triplet basis, consisting of |tt) 
(1,0,0,0^, (l/y2)(|n) + UT)) := (0,1,0,0^, Ui) := 
(0,0,1,0)^ and {l/V2)i\n) - lit)) := (0,0,0,1)^. In 
the case of flux qubits, the t and | states correspond to 
clockwise and counterclockwise current respectively. 

In this basis, the undamped Hamiltonian 1l2qb, equa- 
tion (^, of the two qubit system assumes the matrix 
form 

/ e + K f] -Ajy ' 
_ 1 rj —K rj Ae 
^ ^2 r] K-e At] 
\ -Ary Ae A?/ -K 



H 



(5) 



withe = ei+ea, rj = (Ai-FA2)/y2, Arj = (Ai-Aa)/^^ 
and Ae = ei — e2. From now on, we concentrate on the 
case of equal parameter settings, Ai ~ A2 and ei = e2. 

If we now also express the coupling to the dissipative 
environment in this basis, we find in the case coupling to 
two uncorrelated distinct baths that 



tt26 

"■2qb 



( 



\ 



K 



rj 

-K 7] 
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(6) 



with s = Xi+ X2 and As = Xi— X2- Here the bath me- 
diates transitions between the singlet and triplet states, 
the singlet is not a protected subspace. 

In the case of two qubits with equal parameters that 
are coupled to one common bath we obtain the matrix 



TTlh 

"-2qb 



/ ( - s + K rj 

f] —K r; 

7] K-e + s 

V -K 



(7) 



where s — 2X and As = 0. One directly recognizes that 
compared to (|^) in this case thermalization to the singlet 
state is impeded, because is block-diagonal in the 
singlet and triplet subspaces. The singlet and triplet are 
completely decoupled from each other and in the case of 
one common bath the singlet is also completely decoupled 
from the bath and thus protected from dissipative effects. 
Therefore, a system in contact with one common bath 
that is prepared in the singlet state will never experience 
any decoherence effects. The singlet state is a protected 
subspace (DFS) p2| , although a trivial, one-dimensional 



III. EIGENENERGIES AND EIGENSTATES OF 
THE TWO-QUBIT HAMILTONIAN 

We calculate exact analytical eigenvalues and eigenvec- 
tors of the unperturbed two-qubit system Hamiltonian in 
the aforementioned symmetric case of (^, which reads 



H 
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296 - -2 I rj K-e 
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This Hamiltonian is block-diagonal and the largest block, 
the triplet, is three-dimensional, i.e. it can be diagonal- 
ized using Cardano's formula. Details of that calculation 
are given in Ref. |2^. The case of non-identical qubits is 
more easily handled numerically. 

In the following, \E1), \E2), \E3) and |£'4) denote the 
eigenstates of the two-qubit system. The eigenenergies 
of the unperturbed Hamiltonian m) depend on the three 
parameters K, e and rj. Figure displays the eigenen- 
ergies in more detail for typical experimentally accessi- 
ble values. The parameters that are choosen for the pa- 
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FIG. 1: Plot of the eigenenergies of the eigenstates \E2), 
\E3) and \E4). From upper left to lower right: 1) K = rj = Es 
and e is varied, 2) K = 10 ■ Es, rj = Es and e is varied; 
the inset resolves the avoided level crossing due to the finite 
transmission amplitude -q, 3) i] — e — E^ and K is varied, 4) 
K = e = Es and r] is varied. 



FIG. 2: Plot of the amplitude of the different singlet /triplet 
states of which the eigenstates denoted by \E1), \E2), \E3} 
and I-E4) are composed for the four eigenstates. In all plots e 
is varied and K and r; are fixed to Es- 



rameters e, rj and K in figure |l| correspond to what can 
be reached in flux qubits. They typically assume val- 
ues of a few GHz resembling the parameters of known 
single- and two-qubit experiments in Delft and at 
MIT |2^. Therefore, we will use a characteristic energy 
scale Eg, which is typically Eg — 1 GHz. The corre- 
sponding scales are ts = 1 ns, ujs — 2t: ■ 1 GHz and 
Tg = Vsih/ks) = 4.8 • 10-2 K. Panel 1) shows that for 
large values of e, two of the eigenenergies are degenerate 
(namely for ri,K the states \E1) and \EA) equal the 

states (1/V2)(|n) - lit)) and {1/V2){\U) + lit)), hence 
the eigenenergies are degenerate) while near zero energy 
bias (magnetic frustation / = 1/2) all four eigenenergies 
might be distinguished. Note also that therefore at zero 
energy bias the transition frequency 0^14 = —tun has a 
local maximum, which, as will be shown below, can only 
be accessed via non-symmetric driving. 

If K is set to a big positive value corresponding to large 
ferromagnetic couphng (figure |l|, panel 2), K — 10 ■ Eg) 
the Hamiltonian (^ is nearly diagonal and hence the 
eigenstates in good approximation are equal to the sin- 
glet/triplet basis states. In this case \E3) equals the 
triplet state (l/\/2)(|Ti) + UT>), \E2) and \E4) equal |tT) 
and III) respectively for positive values of e. For large 
negative values of e, the two states \E2) and |i?4) be- 
come equal III) and |tT) with a pseudo-spin-flip between 
clockwise and counter-clockwise rotating current at e = 
when going from positive to negative e. In the case of 
large ferromagnetic coupling, the ground state tends to- 
wards the superposition (l/\/2)(|tt) + III)). Panel 2) 
shows that only for e equal to zero, both \E2) = |tt) 
(|£;2) = III) for negative e) and |£;4) = |||) (|^4) = |TT) 



for negative e) have the same energies (which one would 

expect if the — (1/2) Kai^'' a^z^"^ term in the Hamiltonian 

dominates), because if e is increased the eicri*'' (i=l,2) 
terms in the Hamiltonian change the energy. 

For large antiferromagnetic coupling, | — iir| ^ e, A 
the states |t|) and ||t) are favorable. In this limit, 
the ground state tends towards (1/V2)(|t|) + ||t)) and 
the energy splitting between (l/-\/2)(|t|) -I- ||t)) and 
(l/\/2)(|||) — III)) vanishes asymptotically, leaving the 
ground state nearly degenerate. 

From figure |l|, panel 3), one directly recognizes that the 
singlet eigenenergy crosses the triplet spectrum, which is 
a consequence of the fact that the singlet does not inter- 
act with any triplet states. At zero energy bias (mag- 
netic frustration / = 1/2 for a flux qubit), none of the 
eigenstates equal one of the triplet basis states (e.g. as 
observed for a large energy bias e), they are rather non- 
trivial superpositions. This is elucidated further in the 
next paragraph. The inset of panel 2) depicts the level 
anti-crossing between the eigenenergies of the two states 
I £'2) and |£'4) due to quantum tunneling. 

In general, the eigenstates are a superposition of sin- 
glet/triplet states. Figure ^ shows, how singlet /triplet 
states combine into eigenstates for different qubit param- 
eters. The first eigenstate \E1) equals (l/\/2)(|t|)-||T)) 
for all times while the other eigenstates |£'2), |i?3) and 
I £'4) are in general superpositions of the singlet /triplet 
basis states. For large values of |e|, the eigenstates 
approach the singlet/triplet basis states. In particu- 
lar at typical working points, where e « 5 • A |l^ the 
eigenstates already nearly equal the singlet/triplet basis 
states. Hence, although the anticrossing described above 
corresponds to the anticrossing used in 0, E3| to demon- 
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FIG. 3: Plot of the absolute value of the transition frequencies 
u!32, UJ42 and 0^31. In the left column K = r) = 0.2 • Es and e 
is varied. Right column K = Q.2 ■ E^, t = Es and t) is varied. 




FIG. 4: Plot of the transition frequencies 0^21, 1^41 and 

In the left column K = rj — 0.2 • £"3 and e is varied. Right 

column K = 0.2 ■ Es, e = Es and ri is varied. 



strate Schrodinger's cat states, entanglement is prevalent 
away from the degeneracy point. For an experimental 
proof, one still would have to show that one has success- 
fully prepared coherent couplings by spectroscopically 
tracing the energy spectrum. Note, that for clarity, in 
figure ^ the inter-qubit coupling strength K is fixed to 
a rather high value of Eg that also sets the width of the 
anti-crossing, which potentially can be very narrow. 



A. Spectroscopy 

As a first technological step towards demonstrating co- 
herent manipulation of qubits, usually the transition fre- 
quencies between certain energy levels are probed ^ , 
i.e. the energy differences between the levels. Figures 
^ and ^ depict the transition frequencies between the 
four eigenstates. The transition frequencies are defined 
as LOnrn = [En ~ Em)/^ and uJnm = ^^mn- The transi- 
tions between the singlet state \E1) and the triplet states 
are forbidden in the case of one common bath, due to the 
special symmetries of the Hamiltonian (^) , if the system 
is driven collectively through a time dependent energy 
bias ei{t) = e2it). However, in the case of two distinct 
baths the environment can mediate transitions between 
the singlet state and the triplet states. 

Not all transition frequencies have local minima at e = 
0. The frequencies oju and W34 have local maxima at 
zero energy bias e. This can already be inferred from 
figure |l|, panel 1), the energy of the eigenstate \E4) has 
a local minimum at e = 0. Similarily, the substructure 
of W34 can be understood from figure |l|: the frequency 
tJ34 has a local maximum at e = 0, because of the local 
minimum of the eigenenergy of the state \E4). First, if e 
is increased, the level spacing of \EA) and |£'3) decreases. 
Then, for larger values of e the level spacing of \EA) and 
I £'3) increases again. Thus, the structure observed for 
uj^A around e = emerges in figure ^ 



IV. BLOCH-REDFIELD-FORMALISM 

In order to describe decoherence in the weak damping 
limit, we use the Bloch-Redfield- Formalism It pro- 
vides a systematic way of finding a set of coupled master 
equations which describes the dynamics of the reduced 
(i.e. the reservoir coordinates are traced out) density ma- 
trix for a given system in contact with a dissipative envi- 
ronment and has recently been shown to be numerically 
equivalent to the more elaborate path-integral scheme 
prf . The Hamiltonian of our two qubit system in con- 
tact with a dissipative environment, eqs. (Q) and (^), has 
the generic "system-|-bath" form 



2qb 



Hi 



(9) 



where Hs is a bath of harmonic oscillators and 'H.int in- 
herits the coupling to a dissipative environment. In our 
case the effects of driving are not investigated. In Born 
approximation and when the system is only weakly cou- 
pled to the environment, Bloch-Redficld theory provides 
the following set of equations for the reduced density ma- 
trix p describing the dynamics of the system ||2^, 

Pnm{t) = ~i^^nmPnm{i) " ^ RnmMPkiit), (10) 

kl 

where [j„,„ = [En - E„i)/fi, and max \Re{Rnmki)\ < 
min ItiJnTnl must hold. The Redfield relaxation tensor 

n^rn 
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Rnmki comprises the dissipative effects of the couphng 
of the system to the environment. The elements of the 
Redfield relaxation tensor are given through golden rule 
rates Eal 



Rn 



.(+) 



imnk' 



(11) 



A. Two qubits coupled to two distinct baths 

We now evaluate the Golden rule expressions in 
( pi] ) in the case of two qubits, each coupled to a 
distinct harmonic oscillator bath. Here, Hi{t) = 
ex.Tp{iHBt/h)Hj exp{—iHBt/fi) denotes the coupling be- 
tween system and bath in the interaction picture, and 
the bracket denotes thermal average of the bath degrees 
of freedom. Writing down all contributions gives 
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£mnk 
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where (ji%m {i = 1,2) are the matrix elements of a, 
with respect to the eigenbasis of the unperturbed Hamil- 



(-) 

£mnk ' 



tonian (g) and likewise for F 

We assume ohmic spectral densities with a Drude- 
cutoff. This is a realistic assumption i.e. for electromag- 
netic noise |l3j and leads to integrals in the rates which 
are tractable by the residue theorem. The cutoff fre- 
quency LUc for the spectral functions of the two qubits 
is typically assumed to be the largest frequency in the 
problem, this is discussed further in section IV E 
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— and J2{ij->) = 



1 



(13) 



The dimensionless parameter a describes the strength of 
the dissipative effects that enter the Hamiltonian via the 
coupling to the environment, described by s and As. In 
order for the Bloch-Redfield formalism, which involves a 
Born approximation in the system-bath coupling, to be 
valid, we have to assume ai/2 «C 1. After tracing out 
over the bath degrees of freedom, the rates read 
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^ [AVi(tj„fe) + AV2(t^„fc)] 
X (coth(/3?ia;„fc/2) - 1) 

+ -\ A2M(w„fc,2) + AiA/(a;„fc,l) 
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with A^ 



"z,em"z,nk^ 
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M{n,i) =V duj 



{coth{f3huj/2)Q + uj), (15) 



here V denotes the principal value. Likewise, 
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dhUirn/^) + 1) 



' imnk 

X (coth(/3?ia;£m/2) + 1 
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The rates T^g^^f, and T^^-^^, might be inserted into (|11|) to 
build the Redfield tensor. Note here that for ojnk and 
respectively, the real part of the rates (which 
is responsible for relaxation and dephasing) is of value 
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To solve the set of differential equations ( |10| ) , it is con- 
venient not to use the superoperator notation but either 
to write p as a vector. In general the Redfield equations 
(|lo| ) without driving are solved by an ansatz of the type 
p{t) = B exp{Ri)B^^ p{0) , where Ri is a diagonal matrix. 
The entries of this diagonal matrix are the eigenvalues of 
the Redfield tensor (O), written in matrix form, includ- 



ing the dominating term iLUnm (cf. equation 10). Here, 
the reduced density matrix p — (pn, . . . , P44)"' is written 
as a vector. The matrix B describes the basis change to 
the eigenbasis of Ri , in which Ri has diagonal form. 



B. Two qubits coupled to one common bath 

For the case of two qubits coupled to one common bath 
we perform the same calculation as in the last section, 
which leads to expressions for the rates analogous to eqs. 
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with A = Ai.n,nk = '^i^ln'^i%+'^i]lm'^i%+'^i^}m'^i%- 

^Sm'^JL' and 

rLL = ^AJ(a.,„) (coth(/3?i^,™/2) -t- 1) -f ^ 



xV J duj 




J(^) 



{coth{phuj /2)uJim + Uj) , 



(18) 
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The difference between tfie rates for the case of two dis- 
tinct baths 0, @ are the two extra terms ^m^z nk 

(2) ^(1) 



and al {^a 



They originate when tracing out the bath 
degrees of freedom. In the case of one common bath 
there is only one spectral function, which we also assume 

2 

to be ohmic J{lo) = {aT%Lo)/{l + ^). For w„fc and 
^im ^ respectively the real part of the rates is of the 



value r 



(+) 

Irnnk 



C. Dynamics of coupled flux qubits with 
dissipation 

The dissipative effects affecting the two qubit system 
lead to decoherence, which manifests itself in two ways: 
The system experiences energy relaxation on a time scale 
T/f = (F/j is the sum of the relaxation rates of the 
four diagonal elements of the reduced density matrix; 
F/j = — Ylin and A„ are the eigenvalues of the ma- 
trix i?n,m,n,m, n,m = 1,...,4.), Called relaxation time, 
into a thermal mixture of the system's energy eigenstates. 
Therefore the diagonal elements of the reduced density 
matrix decay to the value given by the Boltzmann fac- 
tors. The quantum coherent dynamics of the system are 
superimposed on the relaxation and decay on a usually 
shorter time scale t^.. = — 1, . . . , 4; j ^ j and 

r^Tim — ~^^Rn'm,n m) termed dephasing time. Thus de- 
phasing causes the off diagonal terms (coherences) of the 
reduced density matrix to tend towards zero. 

First, we investigate the incoherent relaxation of the 
two qubit system out of an eigenstate. At long times, the 
system is expected to reach thermal equilibrium, p^q = 
{1/Z)e^'^^ . Special cases are T ~ where p^q equals 
the projector on the ground state and T oo where all 
eigenstates are occupied with the same probability, i.e. 
Poq = (l/4)i. Figure || and || illustrate the relaxation 
of the system prepared in one of the four eigenstates for 
temperatures T — and T — 21 ■ Eg respectively. The 
qubit energies if, e and 77 are all set to Eg and a is set 
to Of = 10~^. From figure |^ one recognizes relaxation 
into the eigenstate \E2), the ground state for this set of 
parameters. 

At low temperatures (T = 0), we observe that for 
the case of two distinct uncorrelated baths a system pre- 
pared in one of the four eigenstates always relaxes into 
the ground state. In the case of two qubits coupling to 
one common bath, this is not always the case, as can be 
seen in the upper left panels of figs. ^ and |[ This can be 
explained through our previous observation, that the sin- 
glet is a protected subspace: Neither the free nor, unlike 
in the case of distinct baths, the bath-mediated dynamics 
couple the singlet to the triplet space. Moreover, we can 
observe that relaxation to the ground state happens by 
populating intermediate eigenstates with a lower energy 
than the initial state the system was prepared in at t = 
(cf. figure |^). 

For high temperatures [T ~ 21 ■ Eg) the system ther- 
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FIG. 5: Plot of the occupation probability of the four eigen- 
states \E1) , \E2) , |_B3) and \E4) for initially starting in one 
of the eigenstates \E1) (first row), \E2) (second row) or |j53) 
(third row) at T = K. The left column illustrates the case 
of two qubits coupling to one common bath and the right 
column the case of two qubits coupling to two distinct baths. 
The energies K, e and "q are all fixed to Eg. The characteristic 
timescale ts is ts = 1/ E^. 



malizes into thermal equilibrium, where all eigenstates 
have equal occupation probabilities. Again, in the case 
of one common bath, thermalization of the singlet state 
is impeded and the three eigenstates \E2) , \E?>) and \EA) 
have equal occupation probabilites of 1/3 after the relax- 
ation time. 

If the system is prepared in a superposition of eigen- 
states, e.g. |i53) and \EA) as in Fig. ^, which are not in a 
protected subspace, we observe coherent oscillations be- 
tween the eigenstates that are damped due to dephasing 
and after the decoherence time the occupation proba- 
bility of the eigenstates is given by the Boltzmann fac- 
tors. This behaviour is depicted in figure ^. Here for 
a — 10~^, the cases of T = and T — 2.1 ■ Eg are com- 
pared. When the temperature is low enough, the system 
will relax into the ground state \E2) , as illustrated by the 
right column of figure M. Thus the occupation probabil- 
ity of the state {1/V2)[\E3) + \E4)) goes to zero. Here, 
in the case of zero temperature, the decoherence times 
for the case of one common or two distinct baths are of 
the same order of magnitude. The left column illustrates 
the behaviour when the temperature is increased. At 
T = 2.1- Eg the system relaxes into an equally populated 
state on times much shorter than for T = 0. For low 
temperatures the characteristic time scale for dephasing 
and relaxation is somewhat shorter for the case of one 
common bath (r^^/r^'' « 0.9 for a = 10~^). This can 
be explained by observing the temperature dependence 




FIG. 6: Plot of the occupation probability of the four eigen- 
states \E1) , \E2) , |-B3) and \E4) for initially starting in one 
of the eigenstates \E1) (upper row) or \E2) (lower row) at 
T — 21- Ea. The left column illustrates the case of two qubits 
coupling to one common bath and the right column the case 
of two qubits coupling to two distinct baths. The energies K, 
e and i] are all fixed to Es- The characteristic timescale ts is 
ts = 1/E,. 



of the rates shown in figure ^. Though for the case of 
one common bath two of the dephasing rates are zero at 
T = 0, the remaining rates are always slightly bigger for 
the case of one common bath compared to the case of 
two distinct baths. If the system is prepared in a general 
superposition, here \E3) and {EA), nearly all rates be- 
come important thus compensating the effect of the two 
rates which are approximately zero at zero temperature 
and leading to faster decoherence. 

If a and therefore the strength of the dissipative effects 
is increased from a = 10~^ to a = 10~^, the observed 
coherent motion is significantly damped. Variation of a 
leads to a phase shift of the coherent oscillations, due 
to renormalization of the frequencies p^ . However, in 
our case the efi'ects of ren ormalization are very small, as 
discussed in chapter IV E , and cannot be observed in our 
plots. 



D. Temperature dependence of the rates 

Figure |^ displays the dependence of typical dephasing 
rates and the relaxation rate Tn on temperature. These 
decoherence rates are the inverse decoherence times. The 
rates are of the same magnitude for the case of one com- 
mon bath and two distinct baths. As a notable excep- 
tion, in the case of one common bath the dephasing rates 



r, 



V21 



= r 



V12 



go to zero when the temperature is de- 
creased while all other rates saturate for T ^ 0. This 
phenomenon is explained later on. If the temperature 
is increased from Tg = (h/kB) ■ Vs — 4.8 • 10~^ K the 
increase of the dephasing and relaxation rates follows a 
power law dependence. It is linear in temperature T with 
a slope given by the prefactors of the expression in the 
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FIG. 7: Plot of the occupation probability 

-^(i/v^){|£;3> + |E4>) (^) when starting in the initial state 
(l/\/2)(|ii'3) + I-B4)), which is a superposition of eigenstates 
_E3) and \EA). First row shows the behaviour for two qubits 
coupling to two uncorrelated baths. The lower row shows 
the behaviour for two qubits coupled to one common bath. 
The qubit parameters e, r) and K axe set to E^, a is set to 



The inset resolves the time scale of the coherent 



a = 10" 
oscillations. 



Redfield rates that depends on temperature. At temper- 
ature T K, Q.l ■ Eg, the rates show a sharp increase for 
both cases. This roll-off point is set by the characteris- 
tic energy scale of the problem, which in turn is set by 
the energy bias e, the transmission matrix element rj and 
the coupling strength K . For the choice of parameters 
in figure ^ the characteristic energy scale expressed in 
temperature is T k Q.\ ■ Eg. 

Note that there is also dephasing between the singlet 
and the triplet states. When the system is prepared (by 
application of a suitable interaction) in a coherent su- 
perposition of singlet and triplet states the phase evolves 
coherently. Then two possible decoherence mechanisms 
can destroy phase coherence. Firstly, "flipless" dephas- 
ing processes, where {E) remains unchanged. These fli- 
pless dephasing processes are described by the terms for 
^im,^nk ^ in the rates, eqs. ( p^ ) and (p^. Obviously 
these terms vanish for T — > 0, as the low- frequency com- 
ponent of Ohmic Gaussian noise is strictly thermal. Sec- 
ondly, relaxation due to emission of a boson to the bath is 
also accompagnied by a loss of phase coherence. This pro- 
cess in general has a finite rate at T = 0. This explains 
the T-dependence of the rates in the single-bath case: 
\E1) alone is protected from the environment. As there 
are incoherent transitions between the triplet-eigenstates 
even at T = 0, the relative phase of a coherent oscillation 
betweeen \E\) and any of those is randomized and the 



decoherence rates F, 



are finite even at T = 0. As 
a notable exception, \E2), the lowest-energy state in the 
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FIG. 8: Log-log plot of the temperature dependence of the 
sum of the four relaxation rates and selected dephasing rates. 
Qubit parameters K, e and rj are all set to Es and a = 10~^. 
The upper panel shows the case of one common bath, the 
lower panel the case of two distinct baths. At the character- 
istic temperature of approximately 0.1 ■ Ts the rates increase 
very steeply. 



FIG. 9: Plot of the temperature dependence of the sum of 
the four relaxation rates and selected dephasing rates. Qubit 
parameters e and rj are set to 0, K is set to Es, a — 10~^ cor- 
responding to the choice of parameters used for the UXOR 
operation. The upper panel shows the case of one common 
bath, the lower panel the case of two distinct baths. In the 
case of one common bath the system will experience no dissi- 
pative effects at T = 0. 



triplet subspace, can only be flipped through absorption 
of energy, which imphes that the dephasing rate T^^i ^-l^o 
vanishes at low temperature. The described behaviour 
can be observed in figure ^. 

If the parameters e and rj are tuned to zero, thus K 
being the only non-vanishing parameter in the Hamilto- 
nian, all dephasing and relaxation rates will vanish for 
T = in the case of one common bath. This behaviour 
is depicted in figure |^. It originates from the special sym- 
metries of the Hamiltonian in this case and the fact that 
for this particular two qubit operation the system Hamil- 
tonian and the coupling to the bath are diagonal in the 
same basis. This special case is of crucial importance for 
the quantum gate operation as described in ^ and affects 
the gate quality factors. 



same parity. The imaginary part of the Redfield tensor 
is given by 
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E. Renormalization effects 

Next to causing decoherence, the interaction with the 
bath also renormalizes the qubit frequencies. This is 
mostly due to the fast bath modes and can be under- 
stood analogous to the Franck-Condon effect, the Lamb 
shift, or the adiabatic renormalization |Q. Renormal- 
ization of the oscillation frequencies LOnm is controlled by 
the imaginary part of the Redfield tensor [ p^ 

Note that lixiRnmnm = —^T^Rmnmn duc to the fact that 
the correlators in the Golden Rule expressions have the 



where V denotes the principal value and C^^^^ are pre- 
factors defined, in the case of two distinct baths, ac- 



cording to C^ 



2b 

£rnnk 



J2) (2) 



the case of one common bath C|^„j, 



and in 



iA. Here, 

for simplicity, we assumed ai ^ a2 ~ a and thus 
Ji{uj) — J2{oj) = J{uj)- Evaluation of the integral leads 



to the following expression for F 
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FIG. 10: The left plot depicts the ratio of the renormalization 
effects and the corresponding transition frequencies. Param- 
eters: a = 10~^, T = and LOc/'^a is varied between 10^ and 
10^ for several frequencies {u>i2, i^i4 and 0^23) for the case of 
two baths and in the case of 0^23 also for the case of one com- 
mon bath. The parameters for the right plot are a — 10"'^, 
T — 2.1 ■ Es and Lod^a is varied between 10^ and 10^. The 
inset of the left plot shows a log-log plot of the temperature 
dependence of the renormalization effects. Here q = 10""^ 
and ujc ~ 10^"^. Note that for small temperatures the renor- 
malization effects do not de pend on temperature. (This is 



elucidated further in chapter IVE.) The plots are scaled log- 



arithmically to emphasize the logarithmic divergence of the 
renormalization effects with ujn. 



with ci := (/3?ia;„fc)/(27r), ci := (/3?iajc)/(27r). In the case 
.(-) 



of V\J,nk expression is 
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+ C2) + i^{c2) 

(23) 



with ci := (w^m/3fi-)/(27r). The terms in ( p2| ) and ( |2^ ) 
which are hnear in Wc give no net contribution to the 
imaginary part of the Redfield tensor jl^. To illus- 
trate the size of the renormalization effects, the ratio of 
the renormalization effects to the frequencies which are 



renormalized is depicted in figure 10 

If ci and C2 are large and the digamma functions can be 
approximated by a logarithm the resulting expression for 
the renormalization effects will be independent of tem- 
perature. The temperature dependence of ( p2| ) and ( |2^ ) 
at higher temperatures, where ci and C2 are small and 
the renormalization effects are very weak is shown in 
figure |l^. The rates (22) and ( p^ ) diverge logarithmi- 
cally with ojc in analogy to the well-known ultraviolet- 
divergence of the spin Boson model . When compar- 
ing the upper left (T = 0) and upper right (T = 2.1 • E'^) 
panel, one recognizes that for the first case one com- 
mon bath gives somewhat smaller renormalization ef- 
fects than two distinct baths, while in the second case 
for T ^ 2.1 ■ Eg the renormalization effects deviate only 
slightly (see the behaviour for UJ23) and the renormal- 
ization effects are smaller for the case of two distinct 
baths. The effects of renormalization are always very 
small {\l'ni{Rn^m,n,m)/'-^nm\ bclow 1% for our choice of 



parameters) and are therefore neglected in our calcula- 
tions. However, having calculated ( p2[ ) and (p3|), these 
are easily incorporated in our calculation. The case of 
large renormalization effects is discussed in Ref. 

We only plotted the size of the renormalization effects 
for LUi2, ujii and W34, because in general all values of 
LUnk are of the same magnitude and give similar plots. 
The size of the renormalization effects diverges linearly 
with a, the dimensionless parameter which describes the 
strength of the dissipative effects. 

For flux qubits, the cutoff frequency lUc is given by the 
circuit properties. For a typical first order low-pass LR 
filter in a qubit circuit |jl^, one can insert i? = 50 
(typical impedance of coaxial cables) and L « 1 nH (de- 
pends on the length of the circuit lines) into lulr = R/L 
and gets that lulr ^ 5 ■ 10^" Hz. ujlr is the largest fre- 
quency in the problem (see again |l^, chapter 4.5) and 
> t^LR should hold. Then ujc « 10" Hz (=10^ • Es) 
as cutoff frequency is a reasonable assumption. 



V. GATE QUALITY FACTORS 



In the last section IV , we evaluated the dephasing and 
relaxation rates of the two-qubit system that is affected 
by a dissipative environment. Furthermore, we visual- 
ized the dynamics of the two-qubit system. This does 
not yet allow a full assessment of the performance as 
a quantum logic element. These should perform unitary 
gate operations and based on the rates alone, one can not 
judge how well quantum gate operations might be per- 
formed with the two qubit system. Therefore, to get a 
quantitative measure of how our setup behaves when per- 
forming a quantum logic gate operation, one can evaluate 
gate quality factors . The performance of a two-qubit 
gate is characterized by four quantities: the fidelity, pu- 
rity, quantum degree and entanglement capability. The 
fidelity is defined as 



1 



16 

16 ^ 



(24) 



where the density matrix obtained from attempting a 
quantum gate operation in a hostile environment is 
Pg ~ p{^g)i which is evaluated for all initial conditions 
p(0) = Y^l^ {^ln\- The fidelity is a measure of how well 
a quantum logic operation was performed. Without dis- 
sipation the reduced density matrix fpQ after performing 
the quantum gate operation, applying Ug and the inverse 
Uq would equal p(0). Therefore the fidelity for the ideal 
quantum gate operation should be 1. 
The second quantifier is the purity 



1 

16 ^ 



(25) 



which should be 1 without dissipation and 1/4 in a fully 
mixed state. The purity characterizes the effects of de- 
coherence. 
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The third quantifier, the quantum degree, is defined 
as the maximum overlap of the resuhing density matrix 
after the quantum gate operation with the maximaUy 
entangled states, the Bell-states 



Q 



\Pg 



(26) 



where the Bell-states ^P^g are defined according to 
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For an ideal entangling operation, e.g. the UXOR gate, 
the quantum degree should be 1. The quantum degree 
characterizes nonlocality. It has been shown ||3^ that all 
density operators that have an overlap with a maximally 
entangled state that is larger than the value 0.78 Q 
violate the Clauser-Horne-Shimony-Holt inequality and 
are thus non-local. 

The fourth quantifier, the entanglement capability C, 
is the smallest eigenvalue of the partially transposed 
density matrix for all possible unentangled input states 
(see below). It has been shown |34 to be negative 
for an entangled state. This quantifier should be -0.5, 
e.g. for the ideal UXOR, thus characterizing a maxi- 
mally entangled final state. Two of the gate quality fac- 
tors, namely the fidelity and purity might also be calcu- 
lated for single qubit gates |T^. However, entanglement 
can only be observed in a system of at least two qubits. 
Therefore, the quantum degree and entanglement capa- 
bility cannot be evaluated for single qubit gates. 

To form all possible initial density matrices, needed 
to calculate the gate quality factors, we use the 16 un- 
entangled product states |^'i„), j = 1,...,16 defined 
according to |*a)i 1*6)2, {a,b = 1,...,4), with 
1*1) = li), 1*2) = IT), 1*3) = (1/V2)(|i) + IT)), and 
1*4) = (l/v/2)(|i) + i IT)). They form one possible ba- 
sis set for the superoperator vq with p{tG) — vgp{^) 
0, H. The states are choosen to be unentangled for 
being compatible with the definition of C. 



A. Implementation of two-qubit operations 

1. Controlled phase-shift gate 

To perform the UXOR operation, it is necessary to 
be able to apply the controlled phase-shift operation to- 
gether with arbitrary single-qubit gates. In the computa- 
tional basis (|00), |01), |10), |11)), the controlled phase- 
shift operation is given by 



Ucz{<p) = 



/ 1 

10 

10 

VO e"fi 



(29) 



and for </? = tt, up to a global phase factor, 

Ucz cxp [i^f^P] exp Jcri^^) exp {i^a^a^^ 

(30) 

Note that in ( |30| ) only cr^-operations, which commute 
with the coupling to the bath, are needed. The controlled 
phase-shift operation together with two Hadamard gates 
and a single-qubit phase-shift operation then gives the 
UXOR gate. 



2. UXOR gate 

Due to the fact that the set consisting of the UXOR (or 
CNOT) gate and the o ne-q ubit rotations, is complete for 
quantum computation |35[ , the UXOR gate is a highly 
important two-qubit gate operation. Therefore we fur- 
ther investigate the behaviour of the four gate quality 
factors in this case. The UXOR operation switches the 
second bit, depending on the value of the first bit of a 
two bit system. In the computational basis, this opera- 
tion has the following matrix form 



U 



XOR. 



10 
10 
1 
10 



(31) 



Up to a phase factor, the two-qubit UXOR (or CNOT) 
operation can be realized by a sequence of five single- 
qubit and one two-qubit quantum logic operation. Each 
of these six operations corresponds to an appropri- 
ate Hamiltonian undergoing free unitary time evolution 
exp(— (i/?i)Hopt). The single-qubit operations are han- 
dled with Bloch-Redfield formalism, like the two-qubit 
operations. We assume dc pulses (instantaneous on and 
off switching of the Hamiltonian with zero rise time of 
the signal) or rectangular pulses 
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UxOR 



exp 



X exp 
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V 2 



X exp 
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V2 



(32) 



where t/cz(7r) is given by (pO|). This generic implemen- 
tation has been chosen in order to demonstrate the com- 
parison to other coupling schemes as well as for com- 
putational convenience, it is not necessarily the optimum 
scheme for application under cryogenic conditions, where 
slow rise-time AC pulses are preferred. Table | shows 
the parameters we inserted into the one- and two-qubit 
Hamiltonian to receive the UXOR operation. In our case 
we assumed = Eg. However, there is no restriction in 
the use of other values for ^. For a typical energy scale 
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No. 


operation 


parameters 


[Es] 


time 


[s] 


1 


^^P * 2 1^ V2 )) 


£2 = -C, A2 




7~1 = 


72 
2« 


2 


/ ■ TV 1\ 

exp [i^^z) 




T2 = 


1 

2? 


3 


exp (i^o-la^) 

I-' \ 4 z z J 


i^ = e 


= 


1 

4^ 


4 


exp (if cr^) 


£2=C 


T4 = 


1 

4^ 


5 


exp (ijal) 




T5 = 


1 


6 


-p(-i(^)) 


£2 = -C, A2 




T6 = 


2? 



TABLE I: Parameters of the Hamiltonians which are needed 
to perform the UXOR gate operation, only the non-zero pa- 
rameters are listed; ^ = _Es in our case. 




FIG. 11: Pulse sequence needed to perform the quantum 
UXOR operation. Here the elements of the unperturbed 
single- and two-qubit Hamiltonian needed to perform a cer- 
tain operation undergoing free unitary time evolution are 
shown. The dotted horizontal lines denote ^ = 0, and the 
horizontal lines are spaced by |^| = Es. The durations of 
each pulse are not equal in general Tj ^ Ti, i,j — 1, . . . , 6 (cf. 
table r 



of 1 GHz, the resulting times from table || are in the 
nanosecond range. 

To better visualize the pulse sequence needed to per- 
form the quantum UXOR operation, which was already 
given in table |, figure ^ depicts the values of the el- 
ements of the Hamiltonians. Interestingly enough we 
find that for the only two-qubit operation included in the 
UXOR operation e and rj are zero. Thus K is the only 
non-zero parameter and 'H.2qb assumes diagonal form. For 
fiux qubits, implementing the pulse sequence figure |ll| 
involves negative and positive values tuning the mag- 
netic frustration through the qubit loop below or above 
/ — 1/2. Note that e.g. for realistic models of inductively 
coupled flux qubits, it is very difficult to turn on the in- 
teraction Hamiltonian between the two qubits without 
the individual az terms in the Hamiltonian. However, 
for the pulse sequence given in (^^, we might simply 
perform the third, fourth and fifth operation of (|3|) at 
once using only the Hamiltonian with both the individual 
az terms and the inter-qubit coupling. 

To obtain the final reduced density matrix after per- 
forming the six unitary operations (32), we iteratively 



determine the density matrix after each operation with 
Bloch-Redfield theory and insert the attained resulting 
density matrix as initial density matrix into the next 
operation. This procedure is repeated for all possible 



FIG. 12: Log-log plot of the temperature dependence of the 
deviations of the four gate quantifiers from their ideal values 
after performing the controlled phase-shift (CPHASE) gate 
operation. Here the temperature is varied from ~ to 2 • i?s. 
In all cases a = ai = Q2 = 10"'^. 



unentangled initial states given in the last section. We 
inserted no additional time intervals between the op- 
erations. This is usually needed, if one applies Bloch- 
Redfield formalism, because it is known to violate com- 
plete positivity on short timescales. However, we circum- 
vent this problem in our calculations by dropping the 
memory after each operation, when we iteratively calcu- 
late the reduced density matrix. This procedure may lead 
to small inaccuracies as compared to using QUAPI [ pT[ , 
which however should not affect our main conclusions. 



B. Temperature dependence 

1. Controlled phase-shift gate 

We have analyzed the gate quality factors in the cases 
of a common and of two distinct baths respectively. In 
figure |l^ the temperature dependence of the deviations 
of the four gate quality factors from their ideal values 
are depicted as a log-log plot. At temperatures below 
T = 2.5 • 10-2 K « 0.5 • T,, the purity and fidelity are 
clearly higher for the case of one common bath, but if 
temperature is increased above this characteristic thresh- 
old, fidelity and purity are slightly higher for the case of 
two baths. 

In the case of one common bath the fidelity, purity and 
entanglement capability are approaching their ideal value 
1, when temperature goes to zero. This is related to the 
fact that in the case of one common bath all relaxation 
and dephasing rates vanish during the two-qubit step of 
the controlled phase-shift gate due to the special symme- 
tries of the Hamiltonian, when temperature goes to zero 
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as depicted in figure ||. 

The controlled phase-shift operation creates entangle- 
ment. The creation of entanglement is impeded by deco- 
herence effects that vanish when temperature approaches 
zero. Therefore, the entanglement capability exhibits the 
same behaviour as the fidelity and purity. For zero dissi- 
pation (a = 0) the quantum degree has the value 0.5 but 
the entanglement capability is -0.5 thus characterizing a 
maximum entangled state. The reason is that the Bell- 
states, which are generated by the controlled phase-shift 
gate from the input states, result in a basis that is dif- 
ferent from the used basis, but can be transformed using 
only local transformations. 

Furthermore, for finite dissipation, figure |^ shows that 
also for the case of two distinct baths, there are only 
three non- vanishing rates for T —^ 0. The system, being 
prepared in one of the 16 initial states, might relax into 
one of the eigenstates that is an entangled state. 

We observe the saturation of the deviation for the case 
of two baths and can directly recognize the effects of the 
symmetries of the controlled phase-shift operation. For 
given a, the fidelity and purity cannot be increased any- 
more by lowering the temperature in the case of two dis- 
tinct baths. Interestingly enough, we find that for two 
qubits coupling to one common bath, the situation is 
different for temperatures below 0.5 • Tg. Above a tem- 
perature of Ts = 4.8 • 10~^ K, the decrease of the gate 
quality factors shows a linear dependence on tempera- 
ture for both cases of one common or two distinct heat 
baths before it again saturates at about 10^ K « 10^ • T,. 
Finite decoherence effects in the fidelity, purity and en- 
tanglement capability at T = for the case of two dis- 
tinct baths are resulting from the coupling of the sys- 
tem to the environment of harmonic oscillators, which 
(at T = 0) are all in their ground states and can be ex- 
cited through spontaneous emission. But for the case of 
one common bath, the deviation from the ideal fidelity 
goes to zero, when temperature goes to zero. This is due 
to the special symmetries {K is the only nonvanishing 
parameter in the two-qubit operation) of the Hamilto- 
nian, which rules out spontaneous emission. These sym- 
metries are also reflected in the temperature dependence 
of the rates, figure ^. There, for one common bath, all 
rates vanish for T ^ 0. Note that these rates only de- 
scribe the two-qubit part of the operation. However, the 
single-qubit part behaves similarly because the terms in 
the single-qubit Hamiltonian are als oc az- 



2. UXOR gate 

Different to the last section, we now add two single- 
qubit operations (Hadamard gates) to the controlled 
phase-shift operation that do not commute with the 
coupling to the bath. In figure |l3|, the deviations of 
the gate quality factors from their ideal values are de- 
picted as a log-log plot. Again, at temperatures below 
T = 2.5 • 10-2 Kw 0.5 • Es, the purity and fidehty are 
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FIG. 13: Log-log plot of the temperature dependence of the 
deviations of the four gate quantifiers from their ideal val- 
ues after performing the UXOR gate operation. Here the 
temperature is varied from ~ to 2 • In all cases 

Of = ai = 02 = 10"'^. The dotted line indicates the upper 
bound set by the Clauser-Horne-Shimony-Holt inequality. 



higher for the case of one common bath, but if temper- 
ature is increased above this characteristic threshold, fi- 
delity and purity are higher for the case of two baths. 
Note that we have chosen a rather large a, this value can 
substantially be improved by means of engineering . 

The fidelity and purity are clearly higher for the case 
of one common bath, when temperature is decreased be- 
low 0.5 • Tg. This is related to the fact that in the case 
of one common bath all relaxation and dephasing rates 
vanish during the two-qubit-step of the UXOR, due to 
the special symmetries of the Hamiltonian, when tem- 
perature goes to zero as discussed in the last paragraph. 
However, the quantum degree and the entanglement ca- 
pability tend towards the same value for both the case of 
one common and two distinct baths. This is due to the 
fact that both quantum degree and entanglement capa- 
bility are, different than fidelity and purity, not defined 
as mean values but rather characterize the "best" possi- 
ble case of all given input states. This results in the same 
value for both cases. 

In the recent work by Thorwart and Hanggi ||l^, the 
UXOR gate was investigated for a <Ty ^ ® ay^ coupling 
scheme and one common bath. They find a pronounced 
degradation of the gate performance, in particular, the 
gate quality factors only weakly depend on tempera- 
ture. They set the strength of the dissipative effects to 
a — 1Q~*. Their choice of parameters was e « 10 ■ Eg, 
A « I- Eg and K « 0.5- Eg which is on the same order of 
magnitude as the values given in table |. As can be seen 
in figure |l^, we also observe only a weak decrease of the 
gate quality factors for both the case of one common bath 
and two distinct baths in the same temperature range dis- 



cussed by Thorwart and Hanggi, both for a = 10 ^ and 
a — 10"'' and overall substantially better values. This is 
due to the fact that for ay ay coupling, the Hamil- 
tonian does not commute with the coupling to the bath 
during the two-qubit steps of the UXOR pulse sequence. 

We observe the saturation of the deviation for both the 
case of two baths and one common bath. For given a, 
the fidelity and purity can not be increased anymore by 
lowering the temperature, different from the behaviour 
for the controlled phase-shift gate that was discussed in 
the last section. This is due to the application of the 
Hadamard gate whose Hamiltonian does not commute 
with the coupling to the bath. Above a temperature of 
10~^ K, the decrease of the gate quality factors shows a 
linear dependence on temperature for both cases. Here, 
different from the controlled phase-shift gate, we observe 
finite decoherence effects in all four gate quantifiers also 
at T = 0, both for the case of one common or two dis- 
tinct heat baths. These decoherence effects are resulting 
from the coupling of the system to the environment of 
harmonic oscillators, which (at T — 0) are all in their 
ground states and can be excited through spontaneous 
emission as already described above. 

The dotted line in figure |l^ shows that the temperature 
has to be less than about T = 21-i<^s = lKin order to 
obtain values of the quantum degree being larger than 
Q ~ 0.78. Only then, the Clauser-Horne-Shimony-Holt 
inequality is violated and non-local correlations between 
the qubits occur as described in JTzt . 

C. Dependence on the dissipation strength 

The deviations from the ideal values of the gate quanti- 
fiers possess a linear dependence on a as expected. Gen- 
erally (if no special symmetries of the Hamiltonian are 
present) there are always finite decoherence effects also 
at T = 0. Therefore, we can not improve the gate quality 
factors below a certain saturation value, when lowering 
the temperature ||l^, as was also discussed in the last 
section. By better isolating the system from the environ- 
ment and by carefully engineering the environment jT^ 
one can decrease the strength of the dissipative effects 
characterized by a. In order to obtain the desired value 
of 0.999 99 for J', V and Q a needs to be below 
a = 10-^ at T = 0.21 • E,. 

D. Time resolved UXOR operation 

To investigate the anatomy of the UXOR quantum 
logic operation, we calculated the occupation probabil- 
ities of the singlet/triplet states after each of the six op- 
erations, of which the UXOR consists. This time resolved 
picture of the dynamics of the two-qubit system, when 
performing a gate operation, gives insight into details 
of our implementation of the UXOR operation and the 
dissipative effects occuring during the operation. Thus 
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FIG. 14: Time resolved UXOR operation. The system is 
initially prepared in the state |00). Occupation probabilities 
of the singlet /triplet states are shown after completion of a 
time step Ti {i = 1,...,6). For a — ai = 02 = 10~^ and 
T — 21 ■ Es — 1 K clearly deviations from the ideal case can 
be observed. Qubit parameters are set according to table |. 



we are able to completely characterize the physical pro- 
cess, which maps the input density matrix po to pout in 
an open quantum system p2[ . When the system is pre- 
pared in the state = |00) the UXOR operation ( |3l| ) 
does not alter the initial state and after performing the 
UXOR operation the final state should equal the initial 
state III) = 1 00). This can clearly be observed in figure 
[1^ . During the UXOR operation occupation probabilities 
of the four states change according to the individual op- 
erations given in (^2|). At T = 21 ■ Eg, the case of 2 baths 
differs significantly from the case of one common bath. 
After the third operation (the two-qubit operation; only 
there the distinction between one common or two distinct 
baths makes sense) occupation probabilities arc different 
for both environments resulting in a less ideal result for 
the case of two baths. 

In figure |l^, the resulting state after performing the 
UXOR operation always deviates more from the ideal 
value (for a = 0, i.e. no dissipation) for the case of two 
distinct baths, if all other parameters are fixed and set to 
the same values for the both cases. The state -P|oo> is less 
close to the ideal occupation probability 1 and the other 
singlet/triplet states are also less close to their ideal value 
for the case of two distinct baths. The case of two distinct 
baths also shows bigger deviations from the ideal case 
(a = 0) during the UXOR operation (see figure [l^). But, 
if the system is initially prepared in the state |tT) = |11) 
the case of two distinct baths shows bigger deviations 
from the ideal case during the UXOR operation, while the 
resulting state is closer to the ideal case for two distinct 
baths compared to one common bath. 

In figures |lj and [ij it looks like there would be no 
decoherence effects (or at least much weaker decoherence 
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FIG. 15: Time resolved UXOR operation. The system is 
initially prepared in the state |11). Occupation probabilities 
of the singlet /triplet states are shown after completion of a 
time step Ti (i = 1, . . . , 6). For q = ai = 02 = 10~^ and 
T — 21 ■ Es = 1 K deviations from the ideal case can be 
observed. Qubit parameters are set according to table U. 



FIG. 16: Time resolved Purity for the UXOR operation. The 
value of the purity after each time step Ti {i = 1, . . . , 6) is 
shown. Here a = Qi = 03 = 10"^ and T = 2.1 ■ = 100 mK 
(lower panel) or T = (upper panel). Qubit parameters are 
set according to table ^ 



effects) after performing the (first two) single-qubit op- 
erations. However, not all input states are affected by 
the decoherence effects the same way. And when we re- 
gard all possible input states, there are finite decoherence 
effects. This can be explained with figure Figure 
depicts the time resolved purity when performing the 
UXOR operation. We clearly observe that there are finite 
decoherence effects for the first single-qubit operations in 
( ^ ) as well. The difference between the single-qubit and 
two-qubit operations is the steeper decrease of the purity 
due to stronger decoherence in the case of the two-qubit 
operation. The upper panel in figure |l^ depicts the be- 
haviour of the purity for T — 0. Decoherence due to the 
az terms in the Hamiltonian will vanish for T — > in the 
case of one common bath. 



VI. CONCLUSION 

We presented a full analysis of the dynamics and de- 
coherence properties of two solid state qubits coupled to 
each other via a generic type of Ising coupling and cou- 
pled, moreover, either to a common bath, or two inde- 
pendent baths. 

We calculated the dynamics of the system and eval- 
uated decoherence times. From the temperature de- 
pendence of the decoherence rates (figure ||), we con- 
clude that both types of environments show a similar 
behaviour; however, in the case of one common bath, 
two of the decoherence rates are zero, and the remaining 
ones are slightly larger than in the case of two distinct 
baths. This temperature dependence is also reflected in 
the characteristics of the gate quality factors from quan- 



tum information theory, which are introduced as robust 
measures of the quality of a quantum logic operation. We 
illustrate that the gate quality factors depend linearly on 
a, as expected. The time resolved UXOR operation (fig- 
ures H andjll) again illustrates the difference between 
one common and two distinct baths, and moreover we 
observe that single-qubit decoherence effects cx dur- 
ing the UXOR operation are weak. The time scales of 
the dynamics of the coupled two qubit system are com- 
parable to the time scales, which were already observed 
in experiments and discussed in the literature [l^ . 

The question, whether one common bath or two dis- 
tinct baths arc less destructive regarding quantum coher- 
ence can not be clearly answered. For low enough tem- 
peratures, coupling to one common bath yields better 
results. However, when the temperature is increased two 
distinct baths do better; in both temperature regimes, 
though, the gate quantifiers are only slightly different for 
both cases. 

Compared to the work of Thorwart |l^ , the interaction 
part of our model Hamiltonian possesses symmetries (the 
Hamiltonian of the two-qubit operation and the errors 
commute) that lead to better gate quality factors. Fur- 
thermore, analysis of the symmetries and error sources of 
our model system can lead to improved coupling schemes 
for solid state qubits. Milburn and co-workers on the 
other hand focused on comparison of classical and quan- 
tum mechanical dynamics ]T5|] and estimated the deco- 
herence properties of two coupled two-state systems. 

Governale [ p6| determined the decoherence properties 
of two coupled charge qubits whose Hamiltonian differs 
from (H) by the type of inter-qubit coupling, namely 



15 



® (J^^ coupling. However, introducing the quality 
factors gives a measure to judge how certain qubit de- 
signs perform quantum gate operations. 

As a next step, one should consider driving, to be able 
to observe and discuss Rabi oscillations in systems of 
two coupled qubits. It should be investigated, how the 
decoherence properties are modified, if one adds more 
qubits to the system. 
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